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We consider extensions of Lemaitre-Tolman-Bondi (LTB) spacetimes to the dissipative case. For 
doing that we previously carry out a systematic study on LTB. This study is based on two different 
aspects of LTB. On the one hand, a symmetry property of LTB will be presented. On the other 
hand, the description of LTB in terms of some fundamental scalar functions (structure scalars) 
appearing in the orthogonal splitting of Riemann tensor will be provided. We shall consider as 
"natural" generalizations of LTB (hereafter referred to as GLTB) either those metrics admitting 
some similar kind of symmetry as LTB, or those sharing structure scalars with similar dependence 
on the metric. 
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I. INTRODUCTION 

LTB dust models [lJ-Q are among the oldest and most 
interesting solutions to Einstein equations. They de- 
scribe spherically symmetric distribution of inhomoge- 
neous non-dissipative dust (see 0, Q for a detailed de- 
scription of these spacetimes). 

They have been used as cosmological models (see 0- 
Hoj and references therein), in the study of gravitational 
collapse and the problem of the cosmic censor ship [ll|- 
[ItJ , and in quantum gravity [H, EH • 

A renewed interest in LTB has appeared, in relation 
with recent observations of type la supernovae, indicating 
that the expansion of the universe is accelerating. Indeed, 
even if it is true that there is general consensus to invoke 
dark energy as a source of anti-gravity for understand- 
ing the cosmic acceleration, it is also true that a growing 
number of researchers consider that inhomogeneities can 
account for the observed cosmic acceleration, without in- 
voking dark energy (see [20l - [26| and references therein). 

Now, in spite of all their interest, LTB spacetimes 
present an important limitation, namely: they do not 
admit dissipative fluxes. This is a serious shortcoming 
since it is already an established fact that gravitational 
collapse is a highly dissipative process (see (27l - [29f and 
references therein). This dissipation is required to ac- 
count for the very large (negative) binding energy of the 
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resulting compact object (of the order of — 10 53 erg). 

Indeed, it appears that the only plausible mechanism 
to carry away the bulk of the binding energy of the col- 
lapsing star, leading to a neutron star or black hole is 
neutrino emission [301 ] - 

Dissipation processes are usually treated invoking 
two possible (opposite) approximations: diffusion and 
streaming out. 

In the diffusion approximation, it is assumed that the 
energy flux of radiation (as that of thermal conduction) 
is proportional to the gradient of temperature. In this 
regime, an equation of transport should be assumed in 
order to obtain the temperature distribution for each 
model. 

The diffusion approximation is in general very sensible, 
since it applies whenever the mean free path of particles 
responsible for the propagation of energy is very small as 
compared with the typical length of the object, a circum- 
stance found very often in astrophysical scenarios. 

In fact, for a main sequence star such as the sun, the 
mean free path of photons at the centre, is of the order of 
2 cm. Also, the mean free path of trapped neutrinos in 
compact cores of densities about 10 12 g. cm.~ 3 becomes 
smaller than the size of the stellar core [3l|, HH . 

Furthermore, the observational data collected from su- 
pernovae 1987A indicate that the regime of radiation 
transport prevailing during the emission process, is closer 
to the diffusion approximation than to the streaming out 
limit 33]. 

However in many other circumstances, the mean free 
path of particles transporting energy may be large 
enough so as to justify the free streaming approximation. 
Therefore we shall include simultaneously both limiting 
cases of radiative transport (diffusion and streaming out), 
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allowing us to describe a wide range of situations. 

On the other hand, at cosmological scales, even if 
it is true that cold dark matter is non-collisional and 
strongly dominated by rest-mass, so that pressure and 
heat flux terms (of kinetic nature) are negligible, it could 
be interesting to test the stability of conclusions and re- 
sults based on the assumptions above, with respect to 
small deviations from those assumptions. In this sense, 
GLTB's with arbitrarily small (but non- vanishing) dissi- 
pative fluxes, could be very helpful. 

From all of the above, the motivations to generalize 
LTB spacetimes as to admit dissipative fluxes are clearly 
justified. 

Therefore, it is our goal in this manuscript to consider 
possible generalizations of LTB spacetimes to the dissi- 
pative case. As it should be obvious, such generalizations 
are not unique. As a "qualitative" guide in our endeav- 
our we shall look for GLTB's "as similar" as possible 
to LTB's. More precise definitions of what we mean by 
"qualitative" guide and "as similar as possible", will be 
provided later. 

Thus for example, exact solutions to Einstein equations 
describing dissipative geodesic fluids have been found in 
[3^ - l38j , however all of them are shear-free and we know 
that a distinct property of LTB is it shear, accordingly 
we shall search for shearing GLTB. Shearing dissipative 
geodesic fluids may be found in [39| and , though they 
do not become LTB in the non-dissipative case. 

Our search will be based on two different types of ar- 
guments. On the one hand, symmetry arguments. We 
shall find a symmetry property of LTB spacetimes and 
we shall assume that the corresponding generalizations 
(GLTB) describing dissipative dust, share the same kind 
of symmetry. On the other hand, we shall describe LTB 
in terms of some scalar functions which emerge from the 
orthogonal splitting of the Riemann tensor. We shall as- 
sume that two of these scalar functions share the same 
form (with respect to metric functions) in both, LTB and 
GLTB. 

In the specific case of localized configurations we have 
to assume that our fluid distribution is bounded by a 
spherical surface. In order to avoid thin shells on such 
a boundary surface Darmois [41| conditions should be 
imposed. 

We would like to emphasize that, even though some 
specific examples are exhibited, our main goal in this 
work consists in providing different techniques to ob- 
tain exact solutions representing geodesic radiating fluids 
whose properties are in some respect similar to LTB. 



II. FLUID DISTRIBUTION, KINEMATIC AL 
VARIABLES AND BASIC EQUATIONS 

We consider a spherically symmetric distribution of 
geodesic fluid, which may be bounded by a spherical sur- 
face S, or not. The fluid is assumed to be pure dust 
undergoing dissipation in the form of heat flow (diffusion 



approximation) and outgoing null fluid (streaming out 
limit). 

Choosing comoving coordinates the general metric can 
be written (in the case of bounded configurations such a 
line element applies to the fluid inside £) as 

ds 2 = -dt 2 + B 2 dr 2 + R 2 {d0 2 + sin 2 Ode/) 2 ), (1) 

where B and R are functions of t and r and are assumed 
positive. We number the coordinates x° = i, x 1 = r, 
x 2 = 9 and x 3 = </>. Observe that B is dimensionless, 
whereas R has the same dimension as r. Also observe 
that t is the proper time. 

The energy-momentum T a p (inside £ if the system is 
bounded) is assumed to have the form 

T a p = fiVaVp +q a Vp + V a qp + el a lp, (2) 

where \i is the energy density, q a the heat flux, e the 
radiation density, V a the four-velocity of the fluid and 
l a a null four-vector. These quantities satisfy 

V a V a = -l, V a q a = 0, 

l a V a = -1, l a l a = 0. (3) 

Since we have chosen a comoving coordinate system, 
we have 

V a = A~ 1 5%, q a = qB- l 8%, 

l a = A~ 1 5q + B~ x 8i, (4) 

where q is a function of t and r, q a = q\ a and x a a un it 
four- vector along the radial direction, satisfying 

X a X a = h X a V a =0, x a = B- 1 6?. (5) 

It may be more convenient to write ([2]) in the form 

T a p = py a Vp + q (V a Xl3 + XaVfj) + £XaXp, (6) 

with 

p, = [i + e, 
q = q + e. 



A. Einstein equations 
For (TTJ) and ([5]), Einstein equations 

G a [j = 87rT Q/ 3, (7) 

read: 
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8nT 01 = ~8nqB = -2 ( ^ - |^ ) , (9) 



(8) c a /3 may be also written as 



(16) 



(17) 



Next, the mass function m(t, r) introduced by Misner 
and Sharp [43[ (see also Q) is given by 
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where dots and primes denote derivatives with respect to 
t and r respectively. Observe that if e 7^ 0, the dissipa- 
tive dust behaves as an anisotropic fluid with vanishing 
tangential stresses [Hj]. 



B. Kinematical variables and the mass function 

The expansion O is given by 



e-v„- | + a | , 



(12) 



and for the shear we have (remember that the four- 
acceleration vanishes) 



(13) 



where h af3 = g Qj3 + V a Vp. 

Using ([J) we obtain the non-vanishing components of 

m 

a u = Hs 2 a, a 22 = = -\r?cj, (14) 

o sin o 



with 
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We can dehne the velocity U of the collapsing fluid 
as the variation of the areal radius (R) with respect to 
proper time, i.e. 



U = R. 

Then ([TBI) can be rewritten as 



B 



1 + C/2 _ 2m(*,r) 
i? 
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With the above we can express © as 

1 



Anq = E 



3R 



7 (Q-a)' 



From (fT 



and 



= -An(eU + qE)R 2 



m' =4tt ( A + ) R'R 2 - 



Equation (J23J) may be integrated to obtain 
to = J AttR 2 f/t + R'dr 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



(assuming a regular centre to the distribution, so to(0) 
0). We may partially integrate (|24l) to obtain 



3m 47T 



(25) 



C. The exterior spacetime and junction conditions 

In the case of bounded configurations, we assume that 
outside £ we have the Vaidya spacetime (i.e. we assume 
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all outgoing radiation is massless), or in the dissipation- 
less case the Schwarzschild spacetime, described by 



ds 2 



2M(v) 



dv 2 -2dpdv + p 2 (de 2 + 



sin 



(26) 

where M(v) denotes the total mass (which is constant in 
the Schwarzschild case) and v is the retarded time. 

The matching of the full nonadiabatic sphere (includ- 
ing viscosity) to the Vaidya spacetime, on the surface 
r = rs = constant, was discussed in |45| (for the discus- 
sion of the shear-free case see [46[ and [47j), However 
observe that we are now including a null fluid within the 
star configuration. 

Now, from the continuity of the first differential form 
it follows (see 45] for details), 



R = p(v), 



(28) 



and 



dtA j] / 2M(v) | Jp 
dt J \ p dv 



(29) 



Whereas the continuity of the second differential form 
produces 



m(t,r) = M(v), 



(30) 



dt^dJl- 2 -^! 



(27) 



and 



(iv_ b_rt 

2 \ ~R ~ ~B~R 



-B 



R R 
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B 



' RT 
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R. 



(31) 



where = means that both sides of the equation are eval- 
uated on £ (observe a misprint in eq.(40) in [45j and a 
slight difference in notation). 

Comparing (l3~Tj) with © and dTUl) one obtains 



(? = 0. 



(32) 



Thus the matching of (JI]) and ([26| on £ implies (|30|) and 
Also, we have 



s L 

e = 



4irp 2 



(33) 



where ([32"|) has been used and Ls denotes the total lumi- 
nosity of the sphere as measured on its surface, which is 
given by 



L = L, 



1 



2m 
P 



,dp 

dv 



and where 



Lno — 



dM 
dv 



dm j dv 
~~dt \~dt 



(34) 



(35) 



is the total luminosity measured by an observer at rest 
at infinity. 

The boundary redshift zs is given by 



dv 



1 + z, 



with 



dv 
~dt 



R!_ 
B 



R 



(37) 



Therefore the time of formation of the black hole is given 
by 



^ • it ■ E + r = o. 



(38) 



Also observe than from 



L 



and (|3T|) it follows 



L, 



(E + U) 2 ' 
and from p9|) . ([20]). ([29)) and (f37|) 



dv 



(39) 



(40) 



Finally it is worth noticing that in the diffussion ap- 
proximation (e = 0, q 7^ 0) it follows at once from (|33|) 
that the total luminosity (is) vanishes, even though 
there are non-vanishing dissipative fluxes (q) within the 
sphere. This result is an obvious consequence of the dust 
condition (vanishing hydrodynamic pressure). 



D. 



Weyl tensor 



The Weyl tensor is defined through the Ricmann tensor 

(36) 

RaBui the Ricci tensor R a p and the curvature scalar 1Z, 
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The electric part of Weyl tensor is defined by 

E a /3 = CunPuV^V, 

with the following non-vanishing components 

2 2 
En = g-B £, 



(41) 



(42) 



1 



-R 2 £, 



E33 = E 2 2 sin 



(43) 



where 
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Observe that we may also write E a p as: 



Ea/3 = £{XaX/3 ~ 2 h <*P)- 



(44) 



(45) 



Finally, using ©, {TOj, (JTTJI with Jl8j and flH| we 
obtain 



3 m 
fl3 



47T (/2 — e) — £. 



(46) 



E. Evolution equations for the expansion and the 
shear, Bianchi identities and a constraint equation 
for the Weyl tensor 

For the system under consideration two equations de- 
scribing the evolution of the expansion (Raychaudhury) 
and the shear (Ellis [48j . [49]]). respectively, can be easily 
derived (see [50| for details), they read: 



'i ,+ B +2q (B + R} +2 RB= - 



(50) 



Finally, the following constraint equation for the Weyl 
tensor may be derived from the Bianchi identities (e.g. 
see [28| for details) 

[£ c)]' + (£ + 47re)— = -127rgi?~. (51) 



III. STRUCTURE SCALARS 

As we mentioned in the Introduction, part of our anal- 
ysis will be based on a set of scalar functions which ap- 
pear in a natural way in the orthogonal splitting of the 
Riemann tensor (see [52| for details). 

Thus, let us introduce the tensor s Y„ a and X a p (which 
are elements of that splitting [53[ . [Hi), defined by: 



Y 



a/3 



R, 



•ory/36. 



and 



(52) 



(53) 



where R* af3l& = \"r\tp~i&R a p P and r) ep7 s denotes the Levi- 
Civita tensor. 

Tensors Y a @ and X a p may also be expressed through 
their traces and their trace-free parts, as 

Y a fs = ^Y T h a p + Y TF (x a Xi3 ~ ~h a p), (54) 

X a[j = ~X T h al } + X TF (x a X0 ~ ^h a p). (55) 

Then from (^2 ]) -(|S5 ]l and using ©, 0, (glj and fj45]) . 
we obtain 



Y T f — £ - 47T6, 



(56) 



6+ 3 e + 3 47r(A + £) ' 



1 , 2 

2 ' -ecr = 47re-£. 



(47) 



(48) 



The two independent components of Bianchi identities 
for the system under consideration read (see [5l| for de- 
tails) : 

• ,_ X B ~uR q' 5R' n 
/ i+( M + e )- + 2^ + | + 2|^=0, (49) 



X T = 8?r/*, X TF = -£- 4ne. (57) 
Also, combining (|25|) and (l46l) with (l56l) we may write 
4tt 



1", 



TF 



-87T6- 



i? 3 



Thus the scalar Ytf may be expressed through the 
Weyl tensor and the outgoing null fluid or in terms of 
the null radiation, the density inhomogeneity and the 
dissipative variables. It is worth recalling that a link 
between Ytf and the Tolman mass has been established 
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in [52] ■ Also, it has been shown that in the geodesic case, 
Yt f controls the stability of the shearfree condition [5(| ■ 
Scalars ([56"]) - ([5?]) for the line element ([TJ) are 



t R B 

R ~ B' 



2 — — I - 

B + ~R I R 
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TF 
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~B' 



(59) 
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+ R 2 



(61) 



In terms of these scalar functions, equations ([4"T]) , (|4"8"]1 
and (1511). become 



e + ie 2 



-Y r 



1 , 2 

2 1 -0tr = -K 



TF, 



(4tt/2 + X TF )' 



3R' 



X TF + 47rg(6-cr)S. 



(62) 

(63) 
(64) 



Now, from in the non-dissipative case we have 

(65) 

0. On the other 



3R f 

(47T/X + Xtf)' = —X TF , 

H 



implying that if Xtf 
hand if ji' — then 



(X T f)' 



then // 
3i?' 



producing 



X 



TF 



R 

m 

R 3 ' 



X' 



TF, 



(66) 



(67) 



Since Xtf must be regular everywhere within the sphere, 
we must put f(t) — 0, therefore we have fj,' = <^=> 
Xtf = 0. In other words, in absence of dissipation, 
energy density inhomogeneity is controled by the scalar 
Xtf- This last result is also valid for an anisotropic fluid 

m- 

In the case of non-dissipative dust, it can be shown 
that a — implies Xtf = 0. Indeed, in this particular 
case it follows from (|6"3")l that Ytf = 0, which because of 
([56]) and ([57]) produces X T f = 0- 

We shall next analyze in some detail some properties 
of LTB spacetimes. 



IV. LEMAITRE— TOLMAN— BONDI METRIC 

In order to obtain the general form of LTB, we assume 
a geodesic, dissipationless fluid (q = e = 0). Then we 
find after integration of ([S]), 



B(t,r) 



R' 



(l + K(r))V2' 



(68) 



where K is an arbitrary function of r. 

Feeding back ([6"5]) in (TTJ) we obtain the Lemaitre- 
Tohnan-Bondi metric (LTB): 



ds 2 



-dt 2 



(R 



>Y2 



l + /t(r) 



dr 2 + R 2 (d6 2 + sin 2 Odcf ). (69) 



The metric ([69]) is usually associated with an inhomoge- 
neous dust source (and so we shall do here), however it 
is worth mentioning that the most general source com- 
patible with LTB spacetimes is an anisotropic fluid (Q, 

Next, the "Euclidean" condition at the centre re- 
quires that the perimeter of an infinitesimally small circle 
around the centre be given by 2nl where I is the proper 
radius of the circle, given in turn by 



Bdr. 



(70) 



On the other hand, the perimeter of any circle, as it fol- 
lows from (jTJ) is given by 2nR implying that in the neigh- 
borhood of the centre we must impose the condition: 



or, using 



R' = B, 



r=0 „ 
K = 0. 



Now from dTSl and ([6"5]) we obtain 

• n 2m , . 
R = — + n{r), 



R 



and from ([22]) and ([23 



m = m(r), 



(71) 



(72) 



(73) 



(74) 



m' = \ix[iR 2 R!. (75) 
Also, the Bianchi identity (l49l) in this case reads 



producing 



B R 

»(b +2 r 



h(r) 



= 0, 



BR 2 ' 



(76) 



(77) 
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or, using 



1. Proper matter collineation 



/' 



3fe(r)(l + «(r)) 1 / 2 

OR 3 )' 



(78) 



where h(r) is a function of integration. 

Depending on k there are three possible solutions of 

1. k = (Parabolic case) 



It can be shown [561 ] that in the pure dust case, there 
exists a vector field £ such that for any LTB spacetime 



£(T a p — 0, 



(85) 



where £^ denotes the Lie derivative with respect to the 
vector field £, which can be shown to be of the form 



= (2m)VV, ^ 3 =t-t bb (r). (79) 



2. ft > (Hyperbolic case) 

Tfl Tfi 

R = — (coshry - 1), -^(sinh?7 - rj) = t- t bb (r). 

(80) 



3. k < (Elliptic case) 



TO TO 
Jl= (l-COBt,), ^ 



(ri-amrj) =t-t bb {r). (81) 



Where t bb (r) is an integration function of r giving 
the value of time for which R(t, r) = 0, not to con- 
found with the center of symmetry R(t, 0). 

In order to prescribe an explicit model we have to provide 
the three functions k(t), m(r) and t bb (r). However, since 
(1691) is invariant under transformations of the form r = 
r(f), we only need two functions of r. 
Scalars Yp, Ytf, Xtf and Xt for ([69)) are: 



and 



" t ^ = -2|-|= (82) 



R R' 



Xtf 



Xt — 



R 2 R R' 



R 2 RR' R 2 2RR> 



2RR' R 2 



RR' 



R 2 R 2 RR' 



(83) 



(84) 



A. Symmetry properties of LTB spacetimes 

It is a simple matter to check that besides the Killing 
vectors associated with spherical symmetry, LTB space- 
times admit no further Killing vectors. Also, it is not dif- 
ficult to verify that they do not admit conformal Killing 
vectors either. It should be clear that such nonadmit- 
tance refers to the general form of LTB, not to any spe- 
cific solution belonging to LTB spacetime. 

However as we shall see below, a specific kind of sym- 
metry may be ascribed to LTB spacetimes. 



whose components are are given by 

e° = F(t), 



(86) 

(87) 
(88) 



where F(t) is an arbitrary function. Thus for any LTB 
spacetime we can always find a vector field £ satisfying 
HMD- 
It should be stressed that £ is not a Killing vector field, 
i.e. it defines a proper matter collineation. 



V. GENERALIZING LTB TO THE DISSIPATIVE 

CASE 

From the previous section it should be clear that LTB 
spacetimes explicitly exclude dissipative fluxes (either 
in the diffusion approximation or in the streaming out 
limit). 

We shall now approach the problem of extending LTB 
so as to include dissipative fluxes (such spacetimes will 
be referred to as GLTB). As a necessary condition we 
shall require that all GLTB's become LTB in the limit 
when dissipative fluxes vanish. 

Since we are searching for spacetimes as "close" as pos- 
sible to LTB, we shall consider geodesic, shearing, dust 
with q ^ 0. It is worth recalling that the pure dust 
(without dissipation) condition implies that the fluid is 
geodesic, however this is no longer true for the dissipa- 
tive case. Therefore the geodesic condition here is non- 
redundant. 

Then integration of @, produces now: 



with 



B(t,r) 



l + K(t,r) = 



R! 



(1 + K(t,r)) 



1/2 : 



AjrqRdt + C(r) 



(89) 



(90) 



Since in the non-dissipative case 
it follows 



C(r) = (1 + «(r)) 



should become ([68)1 . 



1/2 



(91) 
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and the Euclidean condition at the center, implies 



VI. THE TRANSPORT EQUATION 



C(r) r =° 1. 



(92) 



Feeding back (l89l) into ([TJ and using (l90|) , we obtain 
for the line element 



ds 2 = -dt 2 



(R> 



[J AirqRdt + C{r)} ' 
R 2 (d6 2 + sin 2 



T dr 2 



(93) 



or, using (li?Tl) 

ds 2 = -dt 2 



7 dr 2 



[J 4nqRdt + (1 + wfr)) 1 /*] 
+ R 2 (d0 2 + sin 2 9dcf) 2 ). (94) 

Scalars Yr, Ytf, Xt and A+f for ([M]) are: 

fl fl' K f R' 3 K \ 



Y TF = -- — 



R R' 1 + K \R' 41 + K 
K 



2(1 + A")' 



(95) 



A' 



1?' 



K 



3 X 



41 + K 



2(1 + if) : 



X, 



2i?i?' 



X 



TF 



R 2 



+ R 2 
RR! 



K 
R 2 



RK 



K' 



RR' R{l + K)'' 



A" 



RK 



RR' R 2 2RR' 2R(1+K) 



(96) 



(97) 



(98) 



In order to proceed further and to obtain specific fami- 
lies of solutions we need to impose additional conditions. 
As mentioned before the criteria to select such conditions 
will be dictated by the requirement that the obtained so- 
lution represents the "closest" possible situation to LTB 
spacetime, including dissipative fluxes. 

Of course, this last requirement is still vague enough 
and allows a great deal of possibilities. In this work we 
shall focus on two possible extensions. On the one hand 
we shall consider extensions based on the characteriza- 
tion of LTB spacetimes in terms of the structure scalars. 
On the other hand we shall propose extensions of LTB 
to the dissipative case, based on the symmetry property 
discussed above. However before considering that, let us 
describe the treatment of the transport equation in the 
pure diffusive case. 



In the diffusion approximation (e = 0, q = q), we shall 
use a transport equation derived from the Miiller-Israel- 
Stewart second order phenomenological theory for dissi- 
pative fluids [57H60|. 

Indeed, it is well known that the Maxwell-Fourier law 
for radiation flux leads to a parabolic equation (diffusion 
equation) which predicts propagation of perturbations 
with infinite speed (see |61J-[64| and references therein). 
This simple fact is at the origin of the pathologies [65[ 
found in the approaches of Eckart (6(| and Landau [67| for 
relativistic dissipative processes. To overcome such diffi- 
culties, various relativistic theories with non-vanishing 
relaxation times have been proposed in the past [57i — 
l60l . l68l [69l ] . The important point is that all these the- 
ories provide a heat transport equation which is not of 
Maxwell- Fourier type but of Cattaneo type [70|, leading 
thereby to a hyperbolic equation for the propagation of 
thermal perturbations. 

The corresponding transport equation for the heat flux 
reads 



ph a PVq p , 1 +q a = -Kh a "(X P +Ta p )~ K T 2 (^T< 

(99) 

where K denotes the thermal conductivity, and T and j3 
denote temperature and relaxation time respectively. In 
the case /3 = we recover the Eckart-Landau equation. 

Observe that, due to the symmetry of the problem, 
equation (|99|) only has one independent component, 
which may be written as 

Pq + q=-§T' -\fiqQ + Pjq- (100) 

In the truncated version of the theory, the last term in 
(|99l) is absent (see for example [tTJ ) , and (I100[) becomes 



Finally observe that if e = we obtain from ([5 

„ _ 9(r) 



(101) 



(102) 



1 B 2 R 2 ' 

where g(r) is an arbitrary function. Then using (|102[) in 
POQI) or (jlOip we obtain (up to the function g(r)) the 
temperature distribution of the fluid in terms of metric 
functions. 



VII. EXTENSIONS OF LTB BASED ON 
STRUCTURE SCALARS 

From an inspection of (|52 ]) . ([55 ]) . ([Ml)-® and flM])-®, 
the following remarks are in order: 

• For any geodesic fluid (dissipative or not), the evo- 
lution of O and a is fully controlled by Ytf and 
Y T 
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• In LTB, scalars Yr and Ytf do not contain k. 

• In GLTB, scalars Yt and Ytf differ from their ex- 
pressions in LTB, by the same term. 

Based on the comments above and on the requirement 
of "maximal" similarity between a LTB and its corre- 
sponding GLTB, let us assume that scalars Yt and Ytf 
share the same expression (in terms of R) in both cases. 
Then it follows at once from (g2) and flSSJ} (or ([96]) ), 



K 



R' 3 K 
~R! ~ 4(1 + K) 



l+K 

Integrating (|103[) we obtain 

R'k 1 / 2 



K 



2(1 + K) 



= 0. 



Ci(r), 



(1+^)3/4 

or using <[90j) and (|9T|) 

(S^i?) 1 / 2 ^' 



Ci(r) 



(103) 



(104) 



(105) 



where G\ (r) is an arbitrary integration function. 
Next, integrating (|104j) we obtain 

4 

K + l = 



Ci(r)2/^ + C 3 (r) 



(106) 



where Ci(r) is another integration function, which may 
be related to C(r) as follows. Taking the t derivative of 
(|106p we have 



K 



8C?(r) 



(107) 



W [^(r)/^r+C 2 (r) 
then, combining ([SO]) with ([TUT) it follows 

C?(r) 



2tt<7 



i?(i?') 2 C?(r)/^+C a (r) 



2 ' 



(108) 



finally, comparing (|105[) and (1108)) we also find 



C 2 (r) = 



2(1- AirqRR' 2 ) 
C(r)+ J AirqRdt' 



(109) 



Thus, the assumption above on Yt and Ytf allows us to 
obtain a GLTB from any given LTB. Indeed, from (11081) . 
taking the time and radial dependence of R from any 
given LTB, we obtain the dissipative flux of the corre- 
sponding GLTB, up to the two functions of r, C\ and 
C*2. The former must satisfy the regularity condition 
Ci(0) = 0, and its vanishing brings back to the starting 
"seed" LTB solution. The remaining physical variables 
follow from the field equations. 

Let us now illustrate the method above, with an ex- 
ample inspired in the parabolic subclass of LTB 
Thus, we choose: 



R{t,r) = f(r){T(r)-t?/\ 



(110) 



where f(r) and T(r) correspond to any specific parabolic 
LTB solution. However, care must be exercised in not 
identifying f(r) with the mass function as it happens in 
LTB. ' 

Using (|110p we may write the integral J jwyi a s: 



dt 
(Rff 



(T(r) - tf^dt 



[f>(r)(T(r)-t) + lf(r)T>(r)Y 



(111) 

introducing the variable T(r) — t = r 3 in (jllll) and inte- 
grating we obtain: 



T 4 dr 



(at 3 + by 



i 



(a 2 r 3 + ab) 3 ^5 



2V3 arctan 



l-2f/¥r 



V3 



ill -In' aT+b 



(112) 



where a = f'(r) and b = |/(r)T'(r). 



Feeding back (ll"l"2"j) into (JTDSJ, using ([TTO"]) . we find: 



27rq 



C 2 (r) 



ar 3 +b 
3v^ 



Ci 2 W/(r)- 



2V3 arctan 



In I aT 3 +b 



C a (r)(c 



b) 



(113) 
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Thus, the function C\ controls the magnitude of the 
dissipation and can be chosen as small as desired in the 
case when perturbations of LTB are required. From the 
above it is evident that these GLTB become LTB when 
dissipative fluxes vanish. Also, it should be clear that if 
the physical properties of the "seed" LTB are reasonable, 
so will be the properties of the corresponding GLTB, at 
least for a sufficiently small value of the function C\(r). 
Finally, in the particular case (e = 0) the temperature 
profile of the model can be obtained from (llOOp (or (|101[) ) 
and dUSJ). 

VIII. EXTENSIONS OF LTB BASED ON 
SYMMETRY PROPERTIES 



Next, from (fTTCj) we obtain £° = F(t), as in LTB, and 
(I114P and (|115l) take the form 

F(t)fi + e^ + 2fxF(t)-2qB^ = 0, (119) 
F(t)q + £Y + q ( F(t) + F(t)| +e^+ ft) = 0. 

This last equation can be rewritten in the form 

e [In (qBe)] ' + F(t) [In (qBF(t))} = 0. (121) 
Multiplying (|12ip by qB, it becomes: 

(qB^Y + (qBF) = 0, (122) 



In this section we shall propose another approach to 
obtain GLTB spacetimes. It consists in assuming that 
any GLTB spacetime shares the same symmetry property 
described in section IV for the corresponding LTB. 

Thus, let us assume now, that our GLTB admits a 
proper matter collineation. 

Then, from the energy momentum tensor for dust with 
dissipation ([5]), we obtain, using (1851) . the three equa- 
tions: 



a partial solution of which may be writen as 

qB^ 1 =-i>(t,r), 
qBF(t)=i>'(t,r). 



(123) 
(124) 



From (|117l) and (I124p we have: 

F(t) K R! 



(114) 



8?r 1 + K R 



(125) 



-£°<RV+/4 -a (Vo + £°§ + + ft) + e ^o = o, 

(115) 

ee + ee'-2q^ + 2e ^°|+^ + ft) =0. (116) 

We shall consider two separated subcases, namely: 

• Dissipation in purely diffusion approximation (q ^ 
Q,e = 0,q = q). 

• Dissipation in the streaming out limit (e ^ 0, q = 
0,q = e). 



Then, using (fTU2l) in (1T21)) with flT2"5l) . and taking 
into account, we get 



87r 5 (r) = RR 1 ' 



K 



(1 + JQ3/2' 



(126) 



which is exactly (| 1 18|) . 

More generally, it can be shown that integrability con- 
ditions for (I119P and (|122l) can always be satisfied. In 
the particular case F(t) = 1, we obtain 



A. Diffusion approximation 

From © and l|89p . we can write: 



In other words the existence of a proper matter 
collineation is consistent with the general form of GLTB 
in the diffusion limit. 

We can integrate (|118l) to obtain 



8irqB = 



k r' 



l + K R 

Combining (JMl), (HHH) and ([TTTj) produces 

K 



8ng(r) = RR 



(l + K) 3 / 2 ' 



(117) 



(118) 



(1 



where c(r) is an arbitrary function. 

As a specific example, let us choose for R(t, r) the form 
for the parabolic subclass of LTB given by (|110l) , then we 
obtain 



11 



1 47rg(r) 
(1 + = ~f(r)f'(r) [f'(r)(T(r) - t) + §/(r)T'(r)] 

and introducing this last expression into (I102p we have 



c(r), 



9(r) 



/ a (r)(T(r)-t)a/3 \ /(r)/'(r) 



•c(r) 



/'(r)(T(r)-f) + -/(r)T'(r) 



(128) 



(129) 



where we have used 

So far the GLTB thus obtained, is defined up to four 
functions of r, which by the invariance with respect to the 
transformation of the radial coordinate reduces to three. 
Two of them can be taken exactly as the corresponding 
in the LTB "seed" model, and the remaining one is ob- 
tained via equations (|122|) - (|124j) . from the assumption 
that £ has the same radial dependence as in the LTB 
case. Observe that due to the junction condition (f52"j) we 

have to impose g(r) — 0. Finally, the temperature profile 
can be easily obtained from (|TU0|) (or ([TUT]) ') and (fT29|) . 



B. Streaming out approximation 

In the case q = and e ^ 0, we obtain from Bianchi 
identities 0SJ and §U\) 



From (pg ]) -([Tg5 ]l we find 

£0 

+ £V + 2/<° + 2^ = 0, (136) 

producing 

e = -4 f2F(r,t) +F(r,t)^ + ^IM.) , (i 37 ) 

with £° = F(r,i). 

Following the requirement of "maximal" similarity be- 
tween a LTB and its corresponding GLTB, let us assume 
£° =F(t), in which case flTTT]) . becomes (|55|) . 

Under that condition, we have from (I133p ~ (jl35p 

F(t)(n+e)+£, 1 (pL + e)' + 2(n+e)F(t)-2eB^ = 0, (138) 



B R\ 



producing 



j( r ) 
BR 2 ' 



(130) 



(131) 



exactly as in the LTB case. However, using ([89 |) - (j9"Tj) we 
get in this case 

_ 3j(r) [jAneRdt+(l + K (r))^] 

where j (r) is a function of integration. 

Next, in this approximation (|1 14[) — (|1 161) take the form 

+ e) + £V + e)' + 2(/x + e)£° - 2eB£ 1 = 0, (133) 



-e ( e° + e°| + e 1 ^ + ft ) - = o, (134) 



^ + eV _ 2e | +2e + =0 . (135) 



F(t)c + £V + e ( F(t) + F(t)| + f + ft 



B ^ B 

eBft = 0,(139) 



F(i)e + £y + 2e(F(t)§+^+ftj =0. (140) 



This last equation can be written as: 

F(t)MeB 2 )i + eiH<Bef)}' = 0, (141) 
or in the form 

F(t)(eB 2 j+±(e(Bef)'=0. (142) 

Also, in this case, we can write (|50|) in the form 

[ln(« 2 )]+^H £ fi 2 )f = 0, (143) 



[e{BR) 2 ] + B(eR 2 Y = 0. 



(144) 



The method to obtain GLTB from any given LTB (in 
the streaming out approximation) may now be sketched 
as follows: 
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• Take any specific "seed" LTB with a given vector 
field £ denning a proper matter collineation (since 
F(t) is arbitrary we put it equal to one). 

• Replace the term (eB 2 ) in (|142l) by its expression 
obtained from (|144[) . 

• Assume that the radial dependence of is the same 
as the corresponding to the LTB "seed" solution. 

• From the form of R corresponding to the "seed" 
LTB, integrate (| 142[) with respect to r to obtain 
e. Any remaining arbitrary function of t may be 
found from junction conditions. 

IX. CONCLUSIONS 

We have presented different alternatives for obtain- 
ing spherically symmetric, geodesic, dissipative dust so- 
lutions to Einstein equations (GLTB). All these alter- 
natives are oriented to produce solutions with specific 
similarities to the non-dissipative case (LTB). For doing 
so we have first carried out a study on LTB spacetimes 
based, on the one hand on some symmetry properties of 
LTB, and on the other on their characterization in terms 
of the structure scalars. 

The generalizations based on structure scalars produce 
a dissipative model for each "seed" LTB, such that dissi- 
pative variables are arbitrary only up to a single function 
of r. 



However, the generalizations based on symmetry prop- 
erties have more degrees of freedom, which may be spec- 
ified depending on the specific problem under considera- 
tion. 

For the purely diffusion case the temperature profile 
may be obtained from the transport equation considered 
in section VI. 

At the present time we foresee two interesting applica- 
tions of GLTB's: One the testing of any result obtained 
from any specific LTB against the presence of small but 
nonvanishing dissipative fluxes, such dissipative pertur- 
bations being represented by exact solutions to the Ein- 
stein equations; and two, the study of relaxational effects 
on different aspects of collapse via the transport equa- 
tion. 

Finally, as already stressed in the Introduction, let us 
recall that our intention here is not to produce specific so- 
lutions with a distinct physical interpretation, but rather 
to pave the way for finding such solutions. 
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